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In this paper, we perform a global constraint on the Ricci dark energy model with both 
the flat case and the non-flat case, using the Markov Chain Monte Carlo (MCMC) method 
and the combined observational data from the cluster X-ray gas mass fraction, Supernovae of 
type la (397), baryon acoustic oscillations, current Cosmic Microwave Background, and the ob- 
servational Hubble function. In the flat model, we obtain the best fit values of the param- 
eters in la, 2a regions: Q m0 = 0.29271^323^ 0.0388 > « = 0.38231™^+™^, Age/Gyr = 
13. 48lo i6-o 21 1 H° ~ 69.09^2 37I3 39- m the non-flat model, the best fit parameters are found 
:„ i_ r. r „m- nrlo .o „ — n qnnq+ °3 67 + 0429 n — n qo^c+o. 0386+0. 0521 _ n n9An +o. 0109+0. 0133 
m la, la regions. iz m0 — U.JUUJ_ .0371-0 .0423 1 Q — u - rf840 -o. 0474-0. 0523 > "* — U.U24u_ . 0130-0. 0153 1 

Age/Gyr = 12.54±g %\t° % , H = 72.89t3 05t3.72- Compared to the constraint results in the ACDM 
model by using the same datasets, it is shown that the current combined datasets prefer the ACDM 
model to the Ricci dark energy model. 
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I. INTRODUCTION 

The observation of the Supernovae of type la (SN la) 0, [|[ provides the evidence that the universe is undergo- 
ing accelerated expansion. Jointing current Cosmic Microwave Background (CMB) anisotropy measurement from 
Wilkinson Microwave Anisotropy Probe (WMAP) @, H[ and the data of the Large Scale Structure (LSS) from SDSS 
S @] j one concludes that there exists an exotic energy component with negative pressure, dubbed dark energy, whose 
density accounts for two-thirds of the total energy density in the universe at present. The simplest but most natural 
candidate of dark energy is the cosmological constant A, with the constant equation of state w\ = — 1. In fact, 
the observational data is mostly consistent with the predictions of the ACDM model However, it suffers the 

so-called fine tuning and cosmic coincidence problems. To avoid these problems, dynamic dark energy models are 



considered as an alternative scenario, such as quintessence [10Hl5l |. phantom [16j . quintom 17J and holographic dark 
energy pj, El etc. 

Although more and more dark energy models have been presented, the nature of dark energy is still a conundrum. 
Under such circumstances, the models which are built according to some fundamental principle are more charming, 
such as holographic dark energy model and agegraphic dark energy model [2(| [2l[ . The former one is on the basis 
of holographic principle, and the latter one is derived from taking the combination between the uncertainty relation 
in quantum mechanics and general relativity into account. In this paper, we focus on the holographic dark energy 
model, which is considered as a dynamic vacuum energy and constructed by considering the holographic principle and 
some features of quantum gravity theory. According to the holographic principle, the number of degrees of freedom 
in a bounded system should be finite and has relations with the area of its boundary. By applying the principle to 
cosmology, one can obtain the upper bound of the entropy contained in the universe. For a system with size L and 
UV cut-off A without decaying into a black hole, it is required that the total energy in a region of size L should not 
exceed the mass of a black hole of the same size, thus L 3 p\ < LM^. The largest L allowed is the one saturating this 
inequality, thus 

3c 2 Ml, 

p^^t 1 ' (i) 

where c is a numerical constant and M p i is the reduced Planck Mass M p i = l/V8nG. It just means a duality between 
UV cut-off and IR cut-off. The UV cut-off is related to the vacuum energy, and IR cut-off is related to the large scale 
of the universe, for example Hubble horizon, event horizon or particle horizon as discussed by [i~8l. Il9j. 
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In the paper [l9|, the author took the future event horizon 

f°° dt' f°° da' 

R M= a J t ^T)= a J a ^ ( 2 ) 

as the IR cut-off L. Although this model is confronted with the causality problem, as pointed out in [l9j], it can reveal 
the dynamic nature of the vacuum energy and provide a solution to the fine tuning and cosmic coincidence problem. 
When c > 1, c = 1 and c < 1, the holographic dark energy behaves like quintessence, cosmological constant and 
phantom respectively. Therefore, in this model, the value of parameter c plays an important role in determining the 
property of holographic dark energy. Then, a model with holographic dark energy proportional to the Ricci scalar 
was proposed by Gao, et. al. in .22], called the Ricci dark energy (RDE). In that paper [22], it has shown that this 
model can avoid the causality problem and naturally solve the coincidence problem of dark energy after Ricci scalar 
is regarded as the IR cut-off L~ 2 : 

I- 2 = R = -6(H + 2H 2 + 4)- (3) 

ar 

An fascinating aspect in the study of dark energy is that the cosmological parameters in a given model can be 
constrained by the increasing observational data. It is found that the best way to constrain is using the combination 
of a thorough observation. Now we give a brief review of the previous works on the combined observational constraints 
of the Ricci dark energy model. In Ref. l23l . we constrained the parameters fi m o and a using 192 SN la data [|| from 
the ESSENCE [HJ and Gold sets [HH^the CMB shift parameter R from three-year WMAP data [H, and the 
BAO parameter A from SDSS [Hj], obtaining the best-fittings: ft m0 = 0.34 ±0.04 and a = 0.38 ±0.03. Subsequently, 
the authors in Ref. [3(| HJ utilized the latest 307 Union SNIa data from the Supernova Cosmology Project (SCP), 
the updated shift parameter 7? from the five-year WMAP data [13], and the independent form of the SDSS BAO 
parameter A [291 ] to constrain the parameters f2 m o and a, whose best-fit values are given by f2 m o = 0-318^q'q24 and 
a = 0.359±S^! in pj. In Ref. [H, the SDSS BAO parameter A is replaced by the measurement of D v {0.35)/D v (0.2) 
from SDSS [34] to investigate the generalized Ricci dark energy model. Then Li, et. al in [35| used the extended 397 
SNIa data from the Union±CFA3 sample [361] . only the values of [r s (zd) / Dy (0.2), r s (zd) / Dy(0.35)] in the measurement 
of BAO [13] and the Maximum likelihood values of [Ia(z*), R(z*), z*] and their inverse covariance matrix in the 



measurement of CMB 32], getting the best constraint results: fl m Q — 0.304, a = 0.363. 

In this paper, we will revisit the RDE model and make a thorough investigation on this model with a completely 
consistent analysis of the combined observations. Comparing with the previous works, we use the more complete 
combinations of the observational datasets from the X-ray gas mass fraction in clusters of galaxies (CBF) [38| . 
397 SN la [111 data, the BAO measurement on the values of [r s (zd) / Dy (0.2), r s (zd)/ Dy (0.35)] and their inverse 
covariance matrix in J37[, the CMB observation [12] on the Maximum likelihood values of [£a (<£*), R(z*), z*] and their 
inverse covariance matrix and the observational Hubble data at fifteen different redshifts, including the three more 
observational data H(z = 0.24) = 79.69 ± 2.32, H(z = 0.34) = 83.8 ± 2.96, and H(z = 0.43) = 86.45 ± 3.27 in [H and 
the observational data [2(| El| . We carry out the global fitting on the RDE model using the Markov Chain Monte 
Carlo (MCMC) method. In addition, in this paper we do not only perform the constraint on the parameters in the flat 
RDE model, but also in the non-flat RDE model, comparing with the constraint results in the standard concordance 
model by using the same combined datasets. 

The paper is organized as follows. In next section, we briefly review the RDE model. In section III, we describe 
the method and data. After we perform the cosmic observation constraint, the results on the determination of the 
cosmological parameters are presented. The last section is the conclusion. 



II. REVIEW THE RICCI DARK ENERGY MODEL 



In this section, we give a brief review on the general formula in the RDE model. With a Friedmann-Robertson- 
Walker (FRW) metric 

ds 2 = -dt 2 + a 2 (t)\-—— + r 2 (d9 2 + sin 2 9d(j) 2 )], (4) 
1 — kr z 

the Einstein field equation can be written as 
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model parameters 
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not-flat RDE 


non-flat ACDM 
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n 97qn+ a0352 +0.0424 
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- 

n nnin+0. 0106+0. 0135 
U.UU1U_ .0105-0.0139 
1 q 7cr+0.51+0.65 
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70 qc+2.69+3.18 
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Xmm 


562.543 


521.669 


539.763 


521.557 



TABLE I: The data fitting results of the parameters with la, 2a regions in RDE model and the ACDM model, where the 
combined observational data from CBF, SN, BAO and CMB and OHD are used. 



where the parameter k denotes the curvature of space k = 1,0, —1 for closed, fiat and open geometries, respectively, 
H is the Hubble function, and p is the energy density of a general piece of matter, and their subscripts m and R 
respectively correspond to matter component, including the cold dark matter p c d m and baryon matter pt, and RDE. 

As suggested by Gao et. al., the energy density of RDE is proportional to the Ricci scalar, thus whose energy 
density is given as 

p R = ZaMliH + 2H 2 + — ) oc R, (6) 

where a is the dimensionless parameter in RDE model, which can be determined though cosmic observation con- 
straints. After changing the variable from the cosmic time t to x = In a, we can rewritten the Friedmann Eq. (5) 
as 

1 1/7 TJ^ 

H " = 3M^ Pm ° e ~ 3X + {a ~ l)ke ~ 2X + a{ 2^x~ + 2R2) - (7) 
pi 

With the help of the definitions as follows: 



E - 77-7 ^r«o - D fjZ ^ > ^ k - "772 ' ( 8 ) 



Ho'™' SM 2 ^ 2 '"" Hi 



the Eq. (7) can be ulteriorly rewritten as 



E 2 = (1 - a)VL k e- 2x + n m0 e- 3x + a{\^- + 2E 2 ). (9) 

2 ax 



Solving this first order differential equation about E 2 , we can obtain 

E 2 =n k e~ 2x + n mQ e- 3x + -^—n mQ e-* x + foe-^-^ 

2 — a 

=n k e- 2x + n m0 e-' 3x + n R (x), (10) 

where / is the integral constant and can be derived by the initial condition E(x = 0) = 1, i.e. fl k + fl m o +&rq = 1 is 
used, which is /_ = 1 — il k — ^z^^mo, ar >d ilu is the definition of the dimensionless RDE density, with the expression 
being 

Ms) = 7^n moe - 3x + (1 - n k - -^-r> m0 ) e -( 4 -f (11) 



III. METHOD AND RESULTS 



In our analysis, we perform a global fitting on determining the cosmological parameters using the Markov Chain 
Monte Carlo (MCMC) method. The MCMC method is based on the publicly available CosmoMC package [_2| and 
the modified CosmoMC package [H, H^, |_4j . including the X-ray cluster gas mass fraction. For our models we 
has modified these packages to add the new parameter a with a prior a G [0.1,0.8]. Besides the parameter a, the 
following basic cosmological parameters (Qi>h 2 ,Q c h 2 ,n k )&re also varying with top-hat priors: the physical baryon 
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FIG. 1: 1-D constraints on individual parameters (a, Age/Gyr,£l m o, Ho) and 2-D contours on these parameters with la, 2a 
errors between each other using the combination of the observational data from CBF, SN, BAO, CMB and OHD in the flat RDE 
model. Dotted lines in the 1-D plots show the mean likelihood of the samples and the solid lines are marginalized probabilities 
for the parameters in the flat RDE model p3 ]. 



density £lbh 2 £ [0.005,0.1], the physical dark matter energy density fi c /i 2 £ [0.01,0.99], and the density of space 
curvature £lk £ [—0.1,0.1]. The cosmological parameter tt m o can be derived from the basic parameters above. We 
use a top- hat prior of the cosmic age i.e.lOGyr < to < 20Gyr and impose a weak Gaussian prior on the physical 



baryon density f^ft 2 = 0.022 ± 0.002 from Big Bang nucleosynthesis 45]. Also, in the data fitting process the seven 
parameters (K, 77, 7, & , a b , Sg, a s ) included in the X-ray gas mass fraction f gas are treated as free parameters. As a 
byproduct the best fitting values of these parameters are obtained. And, these values can also be taken accounted as 
a check of data fitting. 

In our calculations, we have taken the total likelihood function L oc e~ x I 2 to be the products of the separate 
likelihoods of CBF, SN, BAO, CMB and OHD. Then we get X 2 is 

X 2 = XCBF + XSN + X 2 BAO + XCMB + XoHD , ( 12 ) 

where the separate likelihoods of CBF, SN, BAO, CMB and OHD and the current observational datasets used in this 
paper are shown in the Appendix |A"1 

The results on the best values of the cosmological parameters with la, 2a errors in the RDE model and the ACDM 
model are listed in Table |TJ In the Fig. [TJ we show one dimensional probability distribution of each parameter 
and two dimensional plots for parameters between each other in the flat RDE model. The corresponding plots in 
the non-flat RDE model arc presented in Fig. [2l From Figs. [T] and [21 it is seen that the cosmological parameters 
in the two cases of RDE model are well determined in la and 2a regions. Comparing the flat RDE model with 
the non-flat RDE model, we can find the difference of Xmin ^ s obvious. The non-flat RDE model with a smaller 
Xmin ~ 539.763 is favored over the flat case with Xmin = 562.543. However, compared to the ACDM model, it is 
found the ACDM models with smaller values Xmin ~ 521.669 (flat case) and Xmin = 521.557 (non-flat case) are better 
fit to the current combined data than the RDE model. What is more, we obtain the best-fit values of the parameters 
in f gas : K = 0.9665, 77 = 0.2058,7 = 1.0866, b = 0.7073, a b = -0.0540, s = 0.1654, a s = 0.1591 in the flat case and 
K = 0.9871, 77 = 0.2114, 7 = 1.0507, b = 0.7749, a b = -0.0950, s = 0.1741, a s = 0.0194 in the non-flat case. 

From Fig. [31 we can see the constraints on the parameters (a, fl m o) m the flat RDE model by using the independent 
dataset from CBF (black solid), SN (red), BAO (blue), CMB (magenta) and OHD (green) and the combined datesets 
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FIG. 2: 1-D constraints on individual parameters (ilk, a , Age/Gyr, H m o, Ho) and 2-D contours on these parameters with la, 2a 
errors between each other using the combination of the observational data from CBF, SN, BAO, CMB and OHD in the non- 
flat RDE model. Dotted lines in the 1-D plots show the mean likelihood of the samples and the solid lines are marginalized 
probabilities for the parameters in the non-flat RDE model [42T ]. 



(black dotted). It is found all the constraint results from the independent dataset are consistent. It is obvious that 
the space of the parameters(a, tt m o) has been reduced by using the three combined datasets from SN, BAO and 
CMB (cyan). It is indicated that the tight constraint has been obtained by the three combined datasets. Combined 
with the additional datasets from CBF and OHD, it is seen that the best fit values of the parameters are changed. 
For the parameters (a, fi TO o), the inclusion of the datasets from CBF and OHD changes their best fit values from 
(0.3630,0.3192) to (0.3823,0.2927). In addition, as shown in Fig. H though the degeneracy with the additional 
datasets from CBF and OHD is not obviously improved, the curves of the 2-D contour plot from the five combined 
datasets become smoother than that from the three combined datasets. 

Next, we investigate the models according to the objective Information Criterion (IC) including the Akaike Infor- 
mation Criterion (AIC) and Bayesian Information Criterion (BIC). Here, we give a brief introduction, for the details 
please see in Ref. [46roC)l |. The AIC was based on the Kullback-Leibler information entropy and derived by H. Akaike. 
It takes the form 

AIC = -2 In £(§\data) max +2K, (13) 

where, C max is the highest likelihood in the model with the best fit parameters 9, K is the number of parameters 
in the model. The first term of Eq. (|13[) measures the goodness of model fit, and the second one interprets model 
complexity. The BIC is similar as AIC, but the second term is different. It was derived by G. Schwarz and is given 
by the form 

BIC = -2 In C(d\data)max + K\nn, (14) 

where, n in the different term is the number of data points in the datasets. In the above two cases, the term 
— 2 In C(9 1 data) max is often called Xmini though it is also generalized to non-Gaussian distributions. 

Now, the problem is how to evaluate which model is the better one. It is the issue of strength of evidence. We take 
AIC case as a example. And, BIC is the same as that. Comparing the AIC values of several models, the minimum 
one is considered as the best value and denoted by AIC m i n = min{AICi,i — 1...N}, where i = 1...N is a set of 
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FIG. 3: The 2-D contours on the parameters (a, f2 m o) with la, 2a and 3a confidence in the flat RDE model for CBF only 
(black solid line), SN only (red line), BAO only (blue line), CMB only (magenta line), OHD only (green line). The joint 
constraint from SN, BAO and CMB is shown in cyan line and the fully combined constraint from CBF, SN, BAO, CMB, and 
OHD is shown in black dotted line. 

alternative candidate models. The difference of AIC for other alternative model i with respect to the besic model is 
expressed as AAICi— AIC i-AIC m i n . The rules of judgement of AIC model selection are that when AAICi is in the 
range — 2, it means the model i has almost the same supports from the data as the best one; for the range 2 — 4, 
this support is considerably less and with AAICi > 10 model i is practically irrelevant [ifjj]. For BIC, the rules read 
that a ABIC of more than 2 (or 6) relative to the best one is considered "unsupported" (or "strongly unsupported") 
from observational data 1471 1 . 



Model 


Parameters 




AIC 


AAIC 


BIC 


ABIC 


flat ACDM 


9 


521.669 


539.669 





576.83 





non-flat ACDM 


10 


521.557 


541.557 


1.888 


582.848 


6.017 


flat RDE 


10 


562.543 


582.543 


42.874 


623.834 


47.003 


non-flat RDE 


11 


539.763 


561.763 


22.094 


607.183 


30.352 



TABLE II: The results of xLn, AIC and BIC in the ACDM model and RDE model and the differences of AIC and BIC with 
respect to the flat ACDM model are listed. 

The values of AIC and BIC in the ACDM model and RDE model and the differences of AIC and BIC with respect 
to the flat ACDM model are listed in Table HI1 In the table, it reads the flat ACDM model are favored by both 
AIC and BIC. According to the rules of judgement of AIC model selection above, it is seen that the non-flat ACDM 
model is supported considerably and the RDE models are disfavored. Since BIC imposes a stricter penalty against 
introducing extra parameters than AIC, as shown from Eq. (fl"4"|) . all the models except the fiat ACDM model are 
disfavored by BIC. 



IV. CONCLUSION 

In summary, in this paper we have performed a global fitting on the cosmological parameters in both the flat 
RDE model and the non-flat RDE model by using a completely consistent analysis of the X-ray gas mass fraction 
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observation, type la supernovae data from the combination of CfA3 samples and the Union set, baryon acoustic 
oscillations data from SDSS, the measurement data on current Cosmic Microwave Background from the five-year 
WMAP observations and the observational Hubble data. The constraint on the parameters in the ACDM model are 
performed by using the same combined datasets. The constraint results are shown in Table U for the flat RDE model 
x 2 mm = 562.543, for the non-flat RDE model xL„ = 539.763, for the flat ACDM model xL„ = 521.669, and for the 
non-flat ACDM model Xmin — 521.557. From Fig. [3l it is shown that the best fit values of the parameters (a, fl m o) 
are changed by the additional datasets from CBF and OHD, though the additional data have a minor effect on the 
confidences of the parameters. It is indicated that more accurate data is anticipated to provide the more stringent 
constraint on the parameters in RDE model. Comparing the RDE model with the ACDM model, we can find that 
the current combined datasets do favor the ACDM model over the RDE model. According to AIC and BIC, we find 
that the flat ACDM model is much supported by the current data. The RDE model is disfavored by AIC and BIC. 
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Appendix A: Cosmological Constraints Methods 

1. The X-ray gas mass fraction constraints 

According to the X-ray cluster gas mass fraction observation, the baryon mass fraction in clusters of galaxies (CBF) 
can be utilized to constrain cosmological parameters. The X-ray gas mass fraction, f gas , is defined as the ratio of 
the X-ray gas mass to the total mass of a cluster, which is approximately independent on the redshift for the hot 
(kT > 5keV), dynamically relaxed clusters at the radii larger than the innermost core f25oo- The X-ray gas mass 
fraction, f gas , can be derived from the observed X-ray surface brightness profile and the deprojected temperature 
profile of X-ray gas under the assumptions of spherical symmetry and hydrostatic equilibrium. Basing on these 
assumptions above, Allen et al. [38[ selected 42 hot (kT > 5keV), X-ray luminous, dynamically relaxed clusters for 
fgas measurements. The stringent restriction to the selected sample can not only reduce maximally the effect of the 
systematic scatter in f gas data, but also ensure that the f gas data is independent on temperature. In the framework 
of the ACDM reference cosmology, the X-ray gas mass fraction is presented as [38[ 



KAjb(z) ( n b 



i + s(z) \n 



d\ CDM (z) 



d A {z) 



1.5 



(Al) 



where A is the angular correction factor, which is caused by the change in angle for the current test model #2500 m 
comparison with that of the reference cosmology O^q 1 ^ 1 - 



A 



_f^ M \\J H(z)d A (z) y 
V feoo J ~ {[H(z)d A (z)} AC ™J ' 



here, the index 77 is the slope of the f ga s(f / '^2500) data within the radius r25oo, with the best-fit average value 
7) = 0.214 ± 0.022 [38|. And the angular diameter distance is given by 

where sinnn(^/|Ofc|a;) respectively denotes sin(-y/|f2fc|a;), v/|f2fe|:c, sinh(y^[fifcja;) for fi^ < 0, fifc = and Qk > 0. 

In equation (jAlj) . the parameter 7 denotes permissible departures from the assumption of hydrostatic equilibrium, 
due to non-thermal pressure support; the bias factor b(z) — b$(\ + &bz) accounts for uncertainties in the cluster 
depletion factor; s(z) = sq(1 + a s z) accounts for uncertainties of the baryonic mass fraction in stars and a Gaussian 
prior for s is employed, with s = (0.16 ± 0.05)/i° 5 [38j]; the factor K is used to describe the combined effects of the 
residual uncertainties, such as the instrumental calibration and certain X-ray modelling issues, and a Gaussian prior 
for the 'calibration' factor is considered by K — 1.0 ± 0.1 (38l |: 
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Following the method in Ref. [38j, l5l| and adopting the updated 42 observational f gas data in Ref. [38| , the best 
fit values of the model parameters for the X-ray gas mass fraction analysis are determined by minimizing, 

N r/ACDMf \ _ f ( \]2 
2 _ Ugas JgasKZi)\ . . , 

XCBF - 2^ ^2-^ . ( A4) 

t /gas v ' 

where 07 oa (^i) is the statistical uncertainties (Table 3 of 38]). As pointed out in _3j| , the acquiescent systematic 
uncertainties have been considered according to the parameters i.e. 77, b(z), s(z) and K . 



2. Type la Supernovae constraints 



We use the 397 SN la Constitution dataset [3g|. The 90 SN la from CfA3 sample with low redshifts are added to 
307 SN la Union sample [52j . The CfA3 sample increases the number of the nearby SN la and reduces the statistical 
uncertainties. In our analysis we use the Constitution datasets used SALT fitter [53| to fit the SN la light curves, 
where the intrinsic uncertainty of 0.138 mag for each CfA3 SNIa, the peculiar velocity uncertainty of 400km/s, and 
the redshift uncertainty of 0.001 have been considered to realize a more cautious assumption that there is the same 
Hubble residual uncertainty between the Cf A3 SN and the nearby Union SN (3(| . Following [54| [55| , one can obtain 
the corresponding constraints by fitting the distance modulus fi(z) as 



Hth(z) = 51og 10 [D L (z)] + m>. 



(A5) 



In this expression Dl(z) is the Hubble-free luminosity distance Hq<Il{z) / c = HocIa(z)(1 + z) 2 /c, with Ho the Hubble 
constant, defined through the re-normalized quantity h as Hq = 100ft. km s~ Mpc~ , and [Xq = 42.38 — 51og 10 h. 



Additionally, the observed distance moduli fj. bs(zi) of SN la at Zi is 

Vobs{zi) = m obs (zi) - M, 



(A6) 



where M is their absolute magnitudes. 

For the SN la dataset, the best fit values of the parameters p s can be determined by a likelihood analysis, based 
on the calculation of 



x 2 {p M') = ^ obs ^ ~ fith ^ Pa ' Zi ^ 2 



SN 



E 

SN 



{51og 10 [D L (p s ,Zi)] - m obs (Zi) +M'Y 



(A7) 



where M' = fio + M is a nuisance parameter which includes the absolute magnitude and the parameter h. The 
nuisance parameter M' can be marginalized over analytically j56|] as 



X 2 {Vs) 



-2 In 



exp 



1 



X 2 (Ps,M') 



dM', 



resulting to 



* B 2 



\u [ — 

2tt 



(A8) 



with 



A = 



B 



C 



E 

SN 

E 

SN 

SN 1 



{5\og 10 [D L (p s ,Zi)} -m obs (zi)} 2 



5log w [D L (p s ,Zi)] -m obs {zi) 
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Relation (|A7[> has a minimum at the nuisance parameter value M' = B/C, which contains information of the values 
of h and M. Therefore, one can extract the values of h and M provided the knowledge of one of them. Finally, note 
that the expression 

X 2 SN {ps,B/C)=A-{B 2 /C), 

which coincides to (|A8[) up to a constant, is often used in the likelihood analysis [13, [H, and thus in this case the 
results will not be affected by a flat M 1 distribution. 



3. Baryon Acoustic Oscillation constraints 

The Baryon Acoustic Oscillations are detected in the clustering of the combined 2dFGRS and SDSS main galaxy 
samples, and measure the distance-redshift relation at z = 0.2. Additionally, Baryon Acoustic Oscillations in the 
clustering of the SDSS luminous red galaxies measure the distance-redshift relation at z = 0.35. The observed scale 
of the BAO calculated from these samples, as well as from the combined sample, are join tly analyzed using estimates 
of the correlated errors to constrain the form of the distance measure Dy{z) [23, [371, [58| 

dz> ^ 1/3 



D ^- c {njm sinAVm L h^) ] ) • (A9) 

The peak positions of the BAO depend on the ratio of Dy(z) to the sound horizon size at the drag epoch (where 
baryons were released from photons) z^, which can be obtained by using a fitting formula [59| : 

1291 CO h 2 ) - 251 

^' ilTmlW»- ll+imhn (A10) 

with 

h = 0.313(fl m /i 2 r°- 419 [l + 0.607(O m /i 2 ) 0674 ], (All) 

b 2 = 0.238(O m /i 2 ) 0223 - (A12) 



In this paper, we use the data of r s (zd)/Dv(z), which are listed in Table llTll where r s (z) is the comoving sound 
horizon size 

* c s dt f a c s da f°° , c. 



r s {z) =c — = c -ttt = c dz 



J a 2 H J z H{z) 

f 1 '^ da 

~VZJo a 2 H{a)yfl + (30 b /(4fi 7 )a) ' 

where c s is the sound speed of the photon— baryon fluid (60l - [62T | : 



(A13) 



-2 n , 9 „ Phi*) _ o . 9 Oft 



4 p-y(z) 4 7 

and here 7 = 2.469 x lQ- b h- 2 for T C mb = 2.725K. 



= 3+-x^=3 + -x(^) a , (A14) 



z 


r s (z d )/Dv(z) 


0.2 


0.1905 ±0.0061 


0.35 


0.1097 ±0.0036 



TABLE III: The observational r s (z d )/D v {z) data ^ 
Using the data of BAO in Table IrU and the inverse covariance matrix V~ x in [37 



V-i = 301241 " 17226 - 9 , (A16) 
\ -17226.9 86976.6 j 
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Thus, the Xbao(Ps) is given as 

X 2 BAO (Ps)=X t V- 1 X, (A16) 
where X is a column vector formed from the values of theory minus the corresponding observational data, with 

/ ^) _ 0.1905 \ 

x = r 1 Z 'V , (A17) 

V 15^5) - 0-1097 ; 

and X 1 denotes its transpose. 

4. Cosmic Microwave Background constraints 

The CMB shift parameter R is provided by [63| 

= ^^sinn[^ f -^-], (A18) 



here, the redshift (the decoupling epoch of photons) is obtained by using the fitting function 64] 

z* = 1048 [l + 0.00124(^ 2 )~ - 738 ] [l + gi (n m h 2 ) 92 ] , 
where the functions g\ and gi read 

, 9l = O.OTSS^^-^^l + SQ.S^ 2 ) - 763 )' 1 , 
.92 = 0.560 (1 + 21. l(^ b /i 2 ) 181 ) _1 . 

In additional, the acoustic scale is related to the distance ratio and is expressed as 



I A = 


7T C 


sinn^]^] / 
Jo 


* dz' 




H( Z >y 




5 — year ML 


5 — year mean 


error, a 


Ia(z*) 


302.10 


302.45 


0.86 


R(z t ) 


1.710 


1.721 


0.019 




1090.04 


1091.13 


0.93 



(A19) 



TABLE IV: The observational l A , R, z* data \3^. 

Using the data of Ia,R,z* in [32j |. which are listed in Table HV| and their covariance matrix of \Ia{z*), R(z*), z*\ 
referring to (321 ]: 



c- 1 



( 1.800 27.968 -1.103 

27.968 5667.577 -92.263 I , (A20) 



\ -1.103 -92.263 2.923 
we can calculate the likelihood L as Xcmb = — 

21nL: 

X 2 CMB = AdilC-^d^d^lAdiY, (A21) 
where Adi = di — df ata is a row vector, and di = (I a, R,z*). 
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5. OHD 

The observational Hubble data are based on differential ages of the galaxies Irlsl. In (6|| , Jimenez et al. obtained an 
independent estimate for the Hubble parameter using the method developed in [651 ] . and used it to constrain the EOS 
of dark energy. The Hubble parameter depending on the differential ages as a function of redshift z can be written in 
the form of 

*<" " -ITlS- (A22 » 

So, once dz/dt is known, H(z) is obtained directly [67|- By using the differential ages of passively-evolving galaxies 
from the Gemini Deep Deep Survey (GDDS) [68[ and archival data [69l - l73 |. Simon et al. obtained H{z) in the range 
of 0.1 < z < 1.8 [13. In SH, Stern et al. used the new data of the differential ages of passively-evolving galaxies at 
0.35 < z < 1 from Keck observations, SPICES survey and VVDS survey. The twelve observational Hubble data from 
US S1I13 are list in Table El Here, we use the value of Hubble constant H = 74.2 ± 3.6km s _1 Mpc _1 , which is 
obtained by observing 240 long-period Cepheids in ■ As pointed out in [13] , the systematic uncertainties have been 
greatly reduced by the unprecedented homogeneity in the periods and metallicity of these Cepheids. For all Cepheids, 
the same instrument and filters are used to reduce the systematic uncertainty related to flux calibration. In addition, 



z 





0.1 


0.17 0.27 


0.4 


0.48 


0.88 


0.9 


1.30 1.43 


1.53 


1.75 


H(z) (km s^Mpc" 1 ) 


74.2 


69 


83 77 


95 


97 


90 


117 


168 177 


140 


202 


la uncertainty 


±3.6 


±12 


±8 ±14 


±17 


±60 


±40 


±23 


±17 ±18 


±14 


±40 



TABLE V: The observational H(z) data [4(J,|4l|]. 

in [3i| , the authors took the B AO scale as a standard ruler in the radial direction, called " Peak Method" , obtaining 
three more additional data: H(z = 0.24) = 79.69 ± 2.32, H(z = 0.34) = 83.8 ± 2.96, and H(z = 0.43) = 86.45 ± 3.27, 
which are model and scale independent. Here, we just consider the statistical errors. 

The best fit values of the model parameters from observational Hubble data (67j are determined by minimizing 

X% U M = £ f g ^^) -f^)] 2 , (A23) 

where p s denotes the parameters contained in the model, H t h is the predicted value for the Hubble parameter, H ^ s 
is the observed value, a(zi) is the standard deviation measurement uncertainty, and the summation is over the 15 
observational Hubble data points at redshifts Zi. 
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